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Abstract 

We present a theorem that allows to simplify linear stability analysis of periodic and quasiperi- 
odic nonlinear regimes in A^-particle mechanical systems (both conservative and dissipative) with 
different kinds of discrete symmetry. This theorem suggests a decomposition of the linearized sys- 
tem arising in the standard stability analysis into a number of subsystems whose dimensions can 
be considerably less than that of the full system. As an example of such simplification, we discuss 
the stability of bushes of modes (invariant manifolds) for the Fermi-Pasta-Ulam chains and prove 
another theorem about the maximal dimension of the above mentioned subsystems. 
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I. INTRODUCTION 



Different dynamical regimes in a mechanical system with discrete-symmetry group Go 
can be classified by subgroups Gj C Gq of this group 

HQS. 

Actually, we can find 

an invariant manifold corresponding to each subgroup Gj and decompose it into the basis 
vectors of the irreducible representations of the group Gq. As a result of this procedure, 
we obtain a bush of modes (see above cited papers) which can be considered as a certain 
physical object in geometrical, as well as dynamical sense. The mode structure of a given 
bush is fully determined by its symmetry group Gj and is independent of the specific type of 
interparticle interactions in the system. In Hamiltonian systems, bushes of modes represent 
dynamical objects in which the energy of initial excitation turns out to be "trapped" (this is 
a phenomenon of energy localization in the modal space). The number of modes belonging 
to a given bush (the bush dimension) does not change in time, while amplitudes of the modes 
do change, and we can find dynamical equations determining their evolution. 

Being an exact nonlinear excitation, in the considered mechanical system, each bush 
possesses its own domain of stability depending on the value of its mode amplitudes. Beyond 
the stability threshold a phenomenon similar to the parametric resonance occurs, the bush 
loses its stability and transforms into another bush of higher dimension. This process is 
accompanied by spontaneous lowering of the bush symmetry: Gj —>■ Gj, where Gj C Gj. 

the detailed theory of these 
Low-dimensional bushes in mechanical systems 

. The 



The concept of bushes of modes was introduced in hJ, 

n 

dynamical objects was developed in 



with various kinds of symmetry and structures were studied in 
problem of bush stability was discussed in 

mm 

. Two last papers are devoted to the 
vibrational bushes in the Fermi-Pasta-Ulam (FPU) chains. 

Note that dynamical objects equivalent to the bushes of modes were recently discussed for 
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the monoatomic chains in the papers of different authors 
that group-theoretical methods developed in our papers 
not only to the monoatomic chains (as was i 
systems with discrete symmetry groups (see. 

In this paper, we present a theorem which can simplify essentially the stability analysis 
of the bushes of modes in complex systems with many degrees of freedom. The usefulness 




3|. Let us emphasize 
can be applied efficiently 

! ), but to all other physical 
)• 



of this theorem is illustrated with the example of nonlinear chains with a large number 
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of particles. Note that the simphfication of the stabihty analysis in such systems actually 
originate from the well-known Wigner theorem about the block-diagonalization of the matrix 
commuting with all matrices of a representation of a given symmetry group. 

In Sec.ini we start with the simplest examples for introducing basic concepts and ideas. In 
Sec. lllH we present a general theorem about invariance of the dynamical equations linearized 
near a given bush with respect to the bush symmetry group. In Sec. II VI we prove a theorem 
which turns out to be very useful for splitting the above mentioned linearized dynamical 
equations for A^-particle monoatomic chains. Some results on the bush stability in the FPU 
chains are discussed in SecfVl 

II. SOME SIMPLE EXAMPLES 

A. FPU-chains and their symmetry 

We consider longitudinal vibrations of A^-particle chains of identical masses (m = 1) and 
identical springs connecting neighboring particles. Let Xi{t) be the displacement of the i- 
th particle {i = 1,2,. ..,N) from its equilibrium position at a given instance t. Dynamical 
equations of such mechanical system (FPU-chain) can be written as follows: 

= f{^i+l ~ ^i) ~ fi^i ~ X'i-l)- (1) 

The nonlinear force f{x) depends on the deformation x of the spring as f{x) = x + x"^ and 
f{x) = X + x^ for the FPU-a and FPU-/5 chains, respectively. We assume the periodic 
boundary conditions 

Xo{t) = XN{t), XN+l{t) = Xi(t) (2) 

to be valid. Let us also introduce the "configuration vector" X{t) which is the A^-dimensional 
vector describing all the displacements of the individual particles at the moment t: 

X{t) = {x,{t),X2{t),...,XN{t)}. (3) 

In the equilibrium state, a given chain is invariant under the action of the operator a 
which shifts the chain by the lattice spacing a. This operator generates the translational 
group 

Tn = {e,a,a'^, . . . ,0,^^'^}, cJ^ = e, (4) 



3 



where e is the identity element and is the order of the cychc group Tj^. The operator 
a induces the cychc permutation of aU particles of the chain and, therefore, it acts on the 
"configuration vector" X{t) as follows: 

aX{t) = a{a;i(t),X2(t),...,Xiv_i(t),XAr(t)} = {xNit), xi{t), X2{t), . . . , XN-i{t)}. 

The full symmetry group of the monoatomic chain contains also the inversion i, with 
respect to the center of the chain, which acts on the vector X{t) in the following manner: 

iX(t) = i{xi(t),X2(t), . . . ,XN-l{t),XN{t)} = {-XAr(t), -XN-l{t), -X2{t), -Xi(t)}. 

The complete set of all products a^i of the pure translations {k = 0,1,2, . . . , N — 1) 
with the inversion i forms the so-called dihedral group Dj^f which can be written as the direct 
sum of two cosets T/v and T/v ■ i: 

= Tjv © Tjv ■ i. (5) 

The dihedral group is a non-Abelian group induced by two generators (a and i) with the 
following generating relations 

= e, i^ = e, la = a~H. (6) 

We will consider different vibrational regimes in the FPU chains, which can be determined 
by the specific forms of the configuration vector. Each of these regimes depends on m 
independent parameters (m < A^) and this number is the dimension of the given regime. 

The simplest case of one- dimensional vibrational regimes represents the so-called vr-mode 
(zone boundary mode) |2l! |: 

X{t) = {A{t), ~A{t) I A{t), -A{t) I A{t), -A{t) I . . . }, (7) 

where A{t) is a certain function of t. In our terminology, this is the one- dimensional bush 
B[a^,i] (see below about notation of bushes of modes). 
The vector 

X{t) = {0, A{t), Bit), 0, -Bit), -Ait) \ ...}, (8) 

represents a two-dimensional vibrational regime that is determined by two time-dependent 
functions Ait) and Bit). This is the two-dimensional bush B[a^,ai] (see j^). 
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In general, for m-dimensional vibrational regime, we write X{t) = C{t), where the 
A^- dimensional vector C{t) depends on m time-dependent functions only. Each specific 
dynamical regime C(t), being an invariant manifold, possesses its own symmetry group that 
is a subgroup of the parent symmetry group Gq = Dj^ of the chain in equilibrium. 



B. FPU-a chain with iV = 4 particles: existence of the bush B[a^,i] 



Let us consider the above discussed equations for the simplest case = 4. Dynamical 
equations ((H) read: 

Xl = f{x2 - Xi) - f{Xi - X4), 
X2 = fix3-X2)- f{x2-Xi), 
X3 = f{Xi - X-i) - f{x3 - X2), 
X4 = f{xi - X4) - f{x4 - X3). 
The symmetry group Gq in the equilibrium state reads: 

Go = Di = [a, i] = {e, d, d^, d'^, i, di, d^i, d'^i}. 

Hereafter, we write generators of any symmetry group in square brackets, while all its 
elements (if it is necessary) are given in curly brackets. 

The operators d and i act on the configuration vector X = {xi,X2,X3,X4} as follows 

aX = {X4, Xi, X2, X3}, iX = {-X4, -X3, -X2, -xi}. 

Therefore, we can associate the following matrices M(d) and M(i) of the mechanical repre- 
sentation with these generators: 

/o 1 \ 
10 
10 
\0 1 0/ 

Their action on the configuration vector X is equivalent to that of the operators d and i, 
respectively. 

Let us now make the transformations of variables in the system (jH)) according to the 
action of the matrices (fTUI) . i.e. 



d M(d) 



i M(i) 














-i\ 










-1 










-1 








V 


-1 








0/ 
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M(d) : Xi 



X4, X2-^Xi, X3-^X2, X4^X3, 



(11) 



M(i) : xi -X4, X2 -X3, X3 -X2, Xi -Xi. (12) 

It is easy to check that both transformations (fTT|l and ()12|1 produce systems of equations 
which are equivalent to the system Q- Moreover, these transformations act on the individual 
equations Uj (j = 1, 2, 3, 4) of the system Q exactly as on the components Xj (j = 1, 2, 3, 4) 
of the configuration vector X. For example, for the operator i (or matrix M(i)) we have 

IMi = — M4, m2 = — M3, lUs = —U2, IU4 = —Ml. 

It is obvious that a certain transposition of these equations multiplied by ±1 will indeed, 
produce a system fully identical to the original system Q- Thus, we are convinced that the 
symmetry group Gq = D4 of our chain in equilibrium turns out to be the symmetry group 
(the group of invariance) of the dynamical equations of this mechanical system. 

Let us now consider the vibrational regime (jZj), i.e. vr-mode, and check that it represents 
an invariant manifold for the dynamical system Q. Substituting xi{t) = x^lt) = A{t), 
X2(t) = X4(t) = —A{t) into (P), we reduce these equations to one and the same equation of 
the form 

A = f{-2A)-fi2A). (13) 

In the case of the FPU-a model this equation turns out to be the equation of harmonic 
oscillator (for the FPU-/5 model it reduces to the Duffing equation) . Indeed, for the FPU-a 
chain, we obtain from the Eq. lfT^ : 

A + AA = 0. (14) 

Using, for simplicity, the initial condition A{0) = Cq, A{0) = 0, we get the following solution 
to Eq. dm): 

A{t) = Co cos{2t). (15) 

Thus, the one-dimensional bush B[a^,i] (or vr-mode) ((Tj) for the FPU-a chain, represents 
purely harmonic dynamical regime 

X{t) = Co{cos(2t), - cos(2t) I cos(2t), - cos(2t)}. (16) 

On the other hand, the invariant manifold X{t) = {A{t), —A{t), A{t), —A(t)}, corre- 
sponding to the bush B[d^, i], can be obtained with the aid of the group-theoretical methods 
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only, without consideration of the dynamical equations (jHI) • Let us discuss this point in more 
detail. 

At an arbitrary instant t, the displacement pattern {A(t), —A(t),A(t), —A{t)} possesses 
its own symmetry group G = D2 G Gq = D4. Indeed, this pattern is conserved under 
inversion (i) and under shifting all particles by 2a. The latter procedure can be considered 
as a result of the action on the chain by the operator a^. These two symmetry elements 
(a^ and i) determine the dihedral group D2 which is a subgroup of order two of the original 
group Go = D, [22]. 

It is obvious, that the old element a of the group Go, describing the chain in equi- 
librium, does not survive in the vibrational state described by the pattern (|7j). Note 
that this element (d) transforms the regime ((71) into its equivalent (but different!) form 



ent 
)■ 



{—A(t),A(t),—A(t),A(t)}. In the present paper, we will not discuss different equiva 
forms of bushes of modes (a detailed consideration of this problem can be found in 
Thus, we encounter the reduction of symmetry Go = -D4 — > G = -D2 when we pass from the 
equilibrium state to the vibrational state ((Zj) for the considered mechanical system. 

The dynamical regime (jZj) represents the one-dimensional bush consisting of only one 
mode (TT-mode). We will denote it as B[G]=B[a^, i] : {A, —A, A, —A}. In square brackets, 
the group of the bush symmetry is indicated by listing its generators (a^ and i, in our case), 
while the characteristic fragment of the bush displacement pattern is presented next to the 
colon. The bush symmetry groiip G fully determines the form (displacement pattern) of the 
bush B[G] (see, for example, 0,0]). Indeed, in the case of the bush B[d^,i], it is easy to show 
that this form, X = {A{t), —A{t),A{t), —A{t)}, can be obtained as the general solution to 
the following linear algebraic equation representing the invariance of the configuration vector 
X: cjiX = X, g2X = X, where gi = and ^2 = 1 are the generators of the group G. In 
our previous papers we often write these invariance conditions for the bush B[G] in the form 

GX = X. (17) 

It is very essential, that the invariant vector X{t), which was found in such geometrical 
(group-theoretical) manner, turns out to be an invariant manifold for the considered dy- 
namical system pf. Thus, we can obtain the symmetry-determined invariant manifolds 
(bushes of modes) without any information on interparticle interactions in the mechanical 
system. 



C. FPU-a chain with = 4 particles: stability of the bush B[d^,i] 



We now turn to the question of the stabihty of the bush B[a^,i], representing a periodic 
vibrational regime X = {A(t), -A{t), A{t), -A{t)}, with A(t) = Cocos(2t). According to 
the conventional prescription, we must linearize the dynamical system Q in the infinitesimal 
vicinity of the given bush and then study the obtained system. For this goal, let us write 



X{t) = C{t)+Sit), 



(18) 



where C = {A{t), —A(t), A(t), —A(t)} represents our bush, while d{t) = 
{Si(t),62(t),63{t),64(t)} is an infinitesimal vector. Substituting (fTSj) into Eqs. (jHl and 
neglecting all terms nonlinear in Sj{t), we obtain the following linearized equations for the 
FPU-a model: 

6i = [62 - 26i + 64] - 4A(t) • [^2 - 54], 
62=[6s- 262 + 6i]+AA{t)-[6s-6i], 
S,= [6,-2d3 + d2]-^A{t)-[6,-62], 
h=[Si-26^ + 6s]+AA{t)-[6i-6s]. 
The last system of equations can be written in the form 



(19) 



d = J(t) • d, 



(20) 



where J(t) is the Jacobi matrix for the system Q calculated by the substitution of the 
vector X = {A{t), —A{t), A{t), —A{t)}. This matrix can be presented as follows: 



J(t) = L + 4A(t) ■ M, 



(21) 



where 



v 



M 



/ 0-1 l\ 
-10 10 
10-1 
\ 1 0-1 / 



(22) 



/-2 1 l\ 
1-210 

1-21 

1 1 -2 / 
are two time-independent symmetric matrices. 

It easy to check that matrices L and M commute with each other: L-M = M-L. Therefore, 
there exists a time-independent orthogonal matrix S that transforms the both matrices L 
and M to the diagonal form: S ■ L ■ S = Ldia, S ■ M ■ S = Mdia (here S is the transposed matrix 
with respect to S). In turn, it means that the Jacobi matrix J(t) can be diagonalized at any 
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S-i 



(23) 



time t by one and the same time- independent matrix S. Therefore, our hnearized system 
(120)) for the considered bush B[a^,i] can be decomposed into four independent differential 
equations. 

Let us discuss how the above matrix S can be obtain with the aid of the theory of 
irreducible representations of the symmetry group G (in our case G = D2). 

In Sec. II VI we will consider a general method for obtaining the matrix S which reduces 
the Jacobi matrix J(t) to a block-diagonal form. This method uses the basis vectors of irre- 
ducible representations of the group G, constructed in the mechanical space of the considered 
dynamical system. In our simplest case of the monoatomic chain with N = 4 particles, this 
method leads to the following result 

/l 1 1 l\ 

1 1-1-1 
1-1 1-1 

VI -1-1 1/ 

The rows of the matrix S from (j^^ are simply the characters of four one-dimensional irre- 
ducible representations (irreps) - Fi, T^, r4 - of the Abelian group D2, because each of 
these irreps is contained once in the decomposition of the mechanical representation of the 
group G = D2- Introducing new variables y = {yi, 2/2, 1/3, 1/4} instead of the old variables 
6 = {6i(t),62(t),6s(t),64^(t)} by the equation y = S ■ 6 with S from we arrive at the 
full splitting of the linearized equations (fTT^ for the FPU-a model: 

m = 0, (24a) 
y2 = -2[l + 4A{t)]y2, (24b) 
2/3 = -%3, (24c) 
y, = -2[l-4A{t)]y,, (24d) 

where A{t) = Cocos(2t). 

With the aid of Eqs. ()24|) . we can find the stability threshold in Cq for loss of stability 
of the one-dimensional bush B[a^,i]. Indeed, according to Eqs. (jSH), the variables yjit) 
{j = 1, 2, 3, 4) are independent from each other, and we can consider them in turn. Eq. ()24a|l 
for yi{t) describes the uniform motion of the center of masses of our chain, since it follows 
from the equations y = S S that yi(t) ~ {Si{t) + 62{t) + 62,{t) + S^it)) . Therefore, considering 
vibrational regimes only, we may assume yi{t) = 0. 



If only y3(t) appears in the solution to the system (j^ . i.e. if yi(t) = 0, y2(t) = 0, y4,(t) = 
0, then we have from the equation 5 = S y (note that S is the orthogonal matrix and, there- 
fore, S"^ = S): d{t) = {y3it),~y3{t),y3it),-y3{t)}, where y^it) ~ cos(2i(;). This solution 
leads only to deviations "along" the bush Xit) = Co{cos(2t), — cos(2t), cos(2t), — cos(2t)} 
and does not signify instability. 

Since A(t) = Cocos(2t), Eq. ()24b|) reads y2 + [2 + 8Co cos(2t)]?/2 = and can be trans- 
formed to the standard form of the Mathieu equation, as well as Eq. ()24d|) . Therefore, the 
stability threshold of the considered bush B[a^,i] for = 4 can be determined directly 
from the well-known diagram of the regions of stable and unstable motion of the Mathieu 
equation. In such a way we can find that critical value Cc for the amplitude Cq of the given 
bush for which it loses its stability is Cc =0.303. 

In conclusion, let us focus on the point that turns out to be very important for proving the 
general theorem in Sec. lIIII The system (fT^ was obtained by linearizing the original system 
(jni), near the dynamical regime C(t) = {A(t), —A{t), A{t), —A{t)}, and Eqs. (jHj) are invariant 
with respect to the parent group Go = [a,i]. Despite this fact, Eqs. (fT!?|l are invariant only 
with respect to its subgroup G = {e, a^, i, a^i} C Go = {e, a, d^, d'^, i, di, d^i, d'^i}: the element 
d G Go (as well as d^, di, d'^i) does not survive as a result of the symmetry reduction Go G. 
Indeed, acting on Eqs. (fT^ by the operator g = a, which transposes variables 6j as follows 

6i 64, 62 61, 63 ^2, h 63, (25) 

we obtain the equations: 

5, = [6i- 26, + 63]-4A{t)-[6,-63], 

6, = [62 - 26i + 54] + AA{t) ■ [62 - 6,], 
= [S, - 262 + 6i]-AA{t)- [63-61], 

6, = [6^- 263 + 62]+AA{t)-[6,-62]. 

Obviously, this system is not equivalent to the system (fnijl ! (The equivalence between (1111)1 
and can be restored, if, besides cyclic permutation in Eqs. (fT^ . we add the artificial 
transformation A(t) — > —A(t)). 

What is the source of this phenomenon? The original nonlinear dynamical system, which 
can be written as X = F{X), is invariant under the action of the operator g = a. Being 
linearized, by the substitution X{t) = C{t) + 5{t) and neglecting all the nonlinear in 6j{t) 
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terms, it becomes 

OF 



dX 



• 5 = J [C{t)] ■ d, (27) 

x=c 

where J [C()f:)] is the Jacobi matrix. The latter system is also invariant under the action of 



the operator g = a, but its transformation must be correctly written as follows 

-1 fdF 
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■ gS = g-'3 [gCit)] ■ gS. (28) 

X=gC 



dX 

In other words, we have to replace the vector X in the Jacobi matrix by a transformed vector, 
gC , near which the linearization is performed. Thus, we must write this matrix in the form 
J [gC{t)] instead of J [C{t)]. In our case, gC{t) = aC{t) = {-A{t),A{t), -A{t),A{t)} and, 
therefore, we indeed have to add the above mentioned artificial transformation A{t) 
-A{t)). 

On the other hand, dealing with the linearized system 6 = J{t) ■ S, we conventionally 
consider the Jacobi matrix J[C(t)] = J(t) as a fixed (but depending on t) matrix which 
does not change when the operator g = a acts on the system S = J{t) ■ S — this operator 
acts on the vector 6 only! The fact is that we try to split the system 6 = J{t) ■ 6 into some 
subsystems using the traditional algebraic transformations of the old variables 6j. Indeed, 
we introduce new variables dnew = S ■ S, where S is a suitable time-independent orthogonal 
matrix, and then obtain the new system Snew = (^S ■ J(t) ■ ■ dnew that decomposes into a 
number of subsystems. 



D. Stability of the bush B[a , i] for the FPU-a chain with > 4 particles 



Linearizing the dynamical equations of the FPU-a chain with = 6 in the vicinity of 
the bush B[d^,i] (vr-mode), we obtain the following Jacobi matrix in Eq. ()20|): 

J(t) = L + AA{t) ■ M, 

where 



/-2 1 l\ 
1-21000 
1-2100 
1-210 
1-21 
\ 1 1 -2 / 



M 



/ -1 l\ 
-10 10 
10-100 
0-1010 
1 -1 
\ 1 0-1 / 
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L-M-M-L 



7^0. 



These two symmetric matrices, unlike the case N = 4, do not commute with each other: 

/ 1 -1 o\ 
-1 1 
-1 1 
10 0-1 
10-1000 
\ -1 1 / 
As a consequence, we cannot diagonahze both matrices L and M simultaneously, i. e. with 
the aid of one and the same orthogonal matrix S. Therefore, it is impossible to diagonalize 
the Jacobi matrix J{t) in equation 6 = J(t) ■ S for all time t. In other words, there are no 
such matrix S that completely splits the linearized system for the bush B[a^,i] for the chain 
with = 6 particles. 

This difference between the cases = 4 and A^ = 6 (generally, for A^ > 4) can be ex- 
plained as follows. The group G = [a^, i] of the considered bush, in fact, determines different 
groups for the cases A^ = 4 and A^ = 6. Indeed, for A^ = 4 [a^,i] = {-E, d^, i, a^i} = D2, 
while for A^ = 6 [d^, i] = {E, a^, d^, i, d^i, d^i} = D^. The latter group (-D3) is non-Abelian 
(id^ = a^i), unlike the group D2 (id^ = a^i) and, as a consequence, it possesses not only 
one-dimensional irreducible representations, but two-dimensional irreps, as well. It will be 
shown in Sec. lIVl that precisely this fact does not permit us to split fully the above discussed 
linearized system 123^ . 

In spite of this difficulty, we can simplify the linearized system d = J{t) ■ 6 considerably 
with the aid of some group-theoretical methods, which are discussed in the two following 
sections. Now, we only would like to present the final result of the above splitting for the 
case A^ = 6: 



yi 



-Ay 



1) 



(29a) 
(29b) 



m + 2i/3 


= Pit)y,, 


2/5 + 22/5 


= my3, 




= Pit)y,, 


y& + 2?/6 


= P{t)y,. 
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(29c) 



(29d) 



Here P{t) = — AA{t)[l + e~^], while P{t) is the complex conjugate function with respect 
to P{t). The two-dimensional subsystems ()29c|) and ()29d|) can be reduces to the real form 
(|87jl (see, Sec. lV A|l by a certain linear transformation. 



Note , th at the stability of the vr-mode ft 



papers 



p, Q Q y, y, Q 



le bush B[a^,i]) was discussed in a number of 
by different methods and with an emphasis 



on different aspects of this stability. In particular, in our paper ^8], a remarkable fact was 
revealed for the FPU-a chain: the stability threshold of the vr-mode is one and the same 
for interactions with all the other modes of the chain. (For other one-dimensional nonlinear 
modes, for both the FPU-a; and FPU-/3 chains, the stability thresholds, determined by 
interactions with different modes, are essentially different 



III. THE GENERAL THEOREM AND ITS CONSEQUENCE 

We consider an A^-degrees-of-freedom mechanical system that described by A^ autonomous 
differential equations 

X = F{X), (30) 

where the configuration vector X = {xi{t),X2(t), . . . ,xj\f(t)} determines the devia- 
tion from the equilibrium state X = {0,0,..., 0}, while vector-function F{X) = 
{fi{X), f2{X), . . . , /n^X)} determines the right-hand-sides of the dynamical equations. 

We assume that Eq. (j30jl is invariant under the action of a discrete symmetry group Gq 
which we call "the parent symmetry group" of our mechanical system. This means that for 
all g E Gq Eq. fl30|) is invariant under the transformation of variables 

X = gX, (31) 

where g is the operator associated with the symmetry element g of the group Gq by the 
conventional definition 

gX = {g~^xi{t), . . . , g-^XNit)}. 
Using dsn and dSH), one can write X = g'^X, g-^X = F{g-'^X), and finally 

X=gF{g-'X). (32) 

On the other hand, renaming X from Eq. (j30| as X, one can write X = F{X). Comparing 
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this equation with Eq. (j32|) . we obtain F{X) = gF{g ^X), or 

FigX) = gFiX). (33) 

This is the condition of invariance of the dynamical equations fl30p under the action of the 
operator g. It must hold for all g E Gq (obviously, it is sufficient to consider such equivalence 
only for the generators of the group Go)- 

Let X{t) = C{t) be an m-dimensional specific dynamical regime in the considered me- 
chanical system that corresponds to the bush B[G] {G C Go). This means that there exist 
some functional relations between the individual displacements Xi{t) (2 = 1,2,..., A^), and, 
as a result, the system ()30p reduces to m ordinary differential equations in terms of the 
independent functions (we denoted them by A{t), B{t), G{t), etc. in the previous section, 
see, for example, Eqs. ()7|8|) ). 

The vector C(t) is a general solution to the equation (see, Eq. (fTTjl ) 

GX = X, 

where G is the symmetry group of the given bush B[G] {G C Gq). 

Now, we want to study the stability of the dynamical regime C{t), corresponding to the 
bush B[G']. To this end, we must linearize the dynamical equations ()3()|1 in a vicinity of the 
given bush, or more precisely, in a vicinity of the vector C{t). Let 

X = C(t) + 5(t), (34) 

where d{t) = {6i{t), . . . , 6N{t)} is an infinitesimal A^-dimensional vector. Substituting X[t) 
from into (jHn|) and linearizing these equations with respect to S{t), we obtain 

5 = J[C{t)] ■ d, (35) 

where J[C(t)] is the Jacobi matrix of the system (IHUj): 



J[C(t)] = 



dfi 




dxj 


X=C{t) 



Now, we intend to prove the following 

Theorem 1. The matrix J[C{t)] of the linearized dynamical equations near a given bush 
B[G], determined by the configuration vector C{t) , commutes with all matrices M{g) (g E G) 
of the mechanical representation of the symmetry group G of the considered bush: 

M{g)-J[C{t)] = J[C{t)]-M{g). 
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Proof. As was already discussed in SeclHl the original nonlinear system X = F{X) trans- 
forms into the system X = g~^F{gX) under the action of the operator g associated with 
the symmetry element g E Gq of the parent group Gq. According to Eq. (jHHj) . the invariance 
of our system with respect to the operator g can be written as follows: 



g-'F{gX) = F{X). 



(36) 



On the other hand, the system X = F{X), linearized in the vicinity of the vector X = C{t) 
reads d = J[C(t)] ■ d (see Eq. flH^ ). Under the action of the operator g, it transforms, 
according to Eq. pH|l . into the system 



5 = g-'j[gC{t)]-gd. 



(37) 



Let us now consider the mechanical representation F of the parent symmetry group Go- 
To this end, we chose the "natural" basis $ = {ci, 62, bat} in the space of all possible 
displacements of individual particles (configuration space): 



ei 








62 



1 





Cat 



/o\ 






(38) 



Acting by an operator g {g E G) on the vector e^, we can write 



N 



j = l,2,...,N. 



(39) 



1=1 



This equation associates the matrix M{g) = ||Mjj|| with the operator g and, therefore, with 
the symmetry element g E G: 

g^g^M{g). (40) 
The set of matrices M{g) corresponding to all g E G forms the mechanical representation F 



for our system 



2J|. As a consequence of this definition, the equation 



gC = M{g) ■ C 



(41) 
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is valid for any vector C determined in the basis (j38p as C = 'Ylik=i ^k^k- 

Using Eq. (jlH), we can rewrite the equation (jHTjl in terms of matrices M{g) = Mg {g G Go) 
of the mechanical representation of the group Gq: 

S = M-^-J [MgC{t)] -Mg-S. (42) 

Therefore, the invariance of the system S = J [C{t)] ■ S with respect of the operator g (matrix 
Mg) can be written as the following relation 

. J [MgCit)] ■Mg = 3 [Cit)] . (43) 

Now, let us suppose that g is an element of the symmetry group G of a given bush B[G] 
[G C Go). By the definition, all the elements of this group {g G G) leave invariant the 
vector C{t) that determines the displacement pattern of this bush 

gCit)=Mg-Cit)=Cit), geG. (44) 

Taking into account this equation, we obtain from the relation 

■ J [C{t)] ■Mg = J [C{t)] , (45) 

which holds for each element g of the symmetry group G of the considered bush. 
Rewriting ()45|) in the form 

J [C{t)] ■ M, = M, ■ J [C{t)] , (46) 

we arrive at the conclusion of our Theorem: all the matrices of the mechanical represen- 
tation of the group G commute with the Jacobi matrix J [C(t)] of the linearized (near the 
given bush) dynamical equations <5 = J [C(t)] • 6. □ 

In what follows, we will introduce a simpler notation for the Jacobi matrix: 

J[C(t)]=J(t). (47) 

Remark. We have proved that all the matrices with g E G commute with the Jacobi 
matrix J(t) of the system ()35|) . But if we take a symmetry element g E Go that is not 
contained in G ((? G Go \ G), the matrix corresponding to g may not commute with J(t). 
An example of such noncommutativity and the source of this phenomenon were presented 
in SecHini 
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Consequence of Theorem^ 



3 



to 



Taking into account Theorem ^ we can apply the well-known Wigner theorem 
split the linearized system 6 = ^{t) ■ 6 into a certain number of independent subsystems. 
Indeed, according to this theorem, the matrix (J(t), in our case) commuting with all the 
matrices of a representation F of the group G (mechanical representation, in our case), can 
be reduced to a very specific block-diagonal form. The dimension of each block of this form 
is equal to Uj * mj, where rij is the dimension of a certain irreducible representation (irrep) 
Tj of the group G containing rrij times in the reducible representation F. Moreover, these 
blocks possess a particular structure which will be considered in Sec. IIVI 

To implement this splitting explicitly one must pass from the old basis ^oid = 
{ei, 62, ... , Bat, } of the mechanical space to the new basis ^new = {<Pi, <p2, ■ ■ ■ , 4>n, } formed 
by the complete set of the basis vectors (pk {k = 1, 2, . . . , iV) of all the irreps of the group 
G. If $„e«, = S ■ ^oid then the unitary transformation 2^ 

Jne«,(t) = S+-J„;d(t)-S (48) 

produces the above discussed block-diagonal matrix ^newif) of the linearized system 5new = 

Jne-!«(^) ■ ^new (here SqIcI S ■ Sfiew^- 

In the next section, we will search the basis vectors ^^^j] (^ = 1;2, ...,nj) of each 
irreducible representation Tj in the form 

04r,] = {x«,xSf,...,xJf)}, (49) 

where x\j\ (A; = 1, 2, ... , N) determines the displacement of the k-th particle corresponding 
to the i-th basis vector ^if^i] of the j-th irrep F^. Actually, this means that we search ^if^j] 
as a superposition of the old basis vectors {k = 1,2, . . . , N) of the mechanical space (see 

my- 

N 

'/>ar,] = ^x(f -e,. (50) 

k=l 

If we find all the basis vectors </>i[rj] in such a form, the coefficients x\j^ are obviously the 
elements of the matrix S that determines the transformation ^new = S ■ ^oid from the old 
basis ^oid = {efc | k = 1, 2, ... , A^} to the new basis ^new = {<Pi[^j] M = 2, . . . , nj; j = 
1,2,...}. Here j = 1, 2, . . . are indices of the irreducible representations that contribute to 
the reducible mechanical representation F. 
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Thus, finding all the basis vectors of the irreps Tj in the form (jl^ provides us 

directly with the matrix S that diagonalizes the Jacobi matrix J(t) of the linearized system 
S = J(t) -(5. 

IV. STABILITY ANALYSIS OF DYNAMICAL REGIMES IN MONOATOMIC 
CHAINS 

A. Setting up the problem and Theorem [2] 

In general, the study of stability of periodic and, especially, quasiperiodic dynamical 
regimes in the mechanical systems with many degrees of freedom presents considerable 
difficulties. Indeed, for this purpose, we must integrate large linearized (near the considered 
regime) system of differential equations with time-dependent coefficients. In the case of 
periodic regime, one can use the Floquet method requiring integration over only one time- 
period to construct the monodromy matrix. But for quasiperiodic regime this method is 
inapplicable, and one often needs to solve system of great number of differential equations 
for very large time intervals to reveal instability (especially, near the stability threshold). 

In such a situation, a decomposition (splitting) of the full linearized system into a number 
of independent subsystems of small dimensions proves to be very useful. Moreover, this de- 
composition can provide valuable information on generalized degrees of freedom responsible 
for the loss of stability of the given dynamical regime for the first time. Let us note that 
the number of such "critical" degrees of freedom can frequently be rather small. 

We want to illustrate the above idea with the case of A^-particle monoatomic chains for 

^ 1. Let us introduce the following notation. The bush B[G] with the symmetry group G 
containing the translational subgroup [a™] will be denoted by B[a™, . . .], where dots stand 
for other generators of the group G. (Note, that any m-dimensional bush can exist only for 
the chain with divisible by m). 

Theorem 2. Linear stability analysis of any bush B[a"^,...] in the N- degrees- of-freedom 
monoatomic chain can be reduced to stability analysis of isolated subsystems of the second 
order differential equations with time- dependent coefficients whose dimensions do not exceed 
the integer number m. 

Corollary. If the bush dimension is d, one can pass on to the subsystems of autonomous 
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differential equations with dimensions not exceeding {m + d). 

Before proving these propositions we must consider the procedure of constructing the 
basis vectors of the irreducible representations of the translational group T. 



B. Basis vectors of irreducible representations of the translational groups 



The basis vectors of irreducible representations o: 



different symmetry groups are usually 
but, in our case, it is easier to make 



obtained by the method of projection operators ^ 

use of the "direct" method based on the definition of the group representation 

Let r be an ra-dimensional representation (reducible or irreducible) of the group G, while 
V[r] be the invariant subspace corresponding to this representation that determined by the 
set $ of A^-dimensional basis vectors (f)j {j = 1, . . . ,n): 

*= {01,02, (51) 

Acting on any basis vector <pj by an operator g {g E G) and bearing in mind the invariance 
of the subspace V[T], we can represent the vector gcpj as a superposition of all basis vectors 
from (jSlj) . In other words, 

^# = {^01,^02,...,^0n} = M((7)#, (52) 

where M{g) is the matrix corresponding, in the representation F, to the element g of the 
group G. (In Eq. (jK^ we use tilde as the symbol of matrix transposition). Eq. (jK^ associates 
with any g E G a. certain n x n matrix M{g) and encapsulates the definition of matrix 
representation 

r = {M((7i),M(^?2),... }. (53) 

The above mentioned "direct" method is based precisely on this definition. Let us use it 
to obtain the basis vectors of the irreducible representations for the translational group T = 
[a™"] . We will construct these vectors in the mechanical space of the A^-particle monoatomic 
chain and, therefore, each vector (pj can be written as follows: 

(f)j = {xi,X2,...,xn}, (54) 
where displacement of i-th. particle from its equilibrium. 
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TABLE I: Irreducible representations of the cyclic group T4 





e 


9 


5^ 


9' 




1 


1 


1 


1 




1 


i 


-1 


—i 




1 


-1 


1 


-1 


r4 


1 


— i 


-1 





The group T . la"] represents a tr^slaU.nal corresponding to the bush 

B[G] =B[a'^, . . .]. For A^-particle chain |22|, T = [a"^] is a subgroup of the order k = N/m 
of the full translational group Tjv = [a] , and we can write the complete set of its elements 
as follows: 

fc = |e, a ,a ,a , . . . , ' \ \a =a = ej. (^ooj 

Being cyclic, the group Tk from (|55|l possesses only one- dimensional irreps, and their total 
number is equal to the order [k = N/m) of this group. 

Below, for simplicity, we consider the case m = 3 and = 12. The generalization to the 
case of arbitrary values of m and turns out to be trivial. 

As it is well-known, the one-dimensional irreps Fj of the /c-order cyclic group can be 
constructed with the aid of fc-degree roots of 1 and, therefore, for our case A^ = 12, m = 3, 
/c = 4, we obtain the irreducible representations listed in Table H] 

In accordance with the definition ()52|1 . the basis vector (p of the one-dimensional irrep F, 
for which M{g) = 7, must satisfy the equation 

g(f) = 70. (56) 

In our case, g = d^, this equation can be written as follows: 

g<p = {Xio, a^ii, Xi2|xi, X2, x^lxi, X5, XqIxt, Xg, Xg} = 

I I I I 1 ^ ^ 

71x1, X2, a;3|x4, X5, xejxr, xg, xg\xio, xu, X12}. 

Here 7 = l,?, —1, —i for the irreps Fi, F2, F3 and F4, respectively. Equating the sequential 
components of both sides of Eq. ()57|1 . we obtain the general solution to the equation gcp = 7^ 



20 



that turns out to depend on three arbitrary constants, say, x, y and z: 

(j) = {x, y, zl'y-^x, 7~^?/, -f~^z\-f~'^x, 7"^?/, 'y-'^z\'y-^x, 7-^?/, 'y~^z} = 
0, 0|7~\ 0, o|7-^ 0, o|7-^ 0, 0}+ 

(58 j 

i/{o,i,o|o, 7-1,010,7-2,010, 7-3,0}+ 

40,0,110,0, 7-^10,0, 7-'|0,0,7-='}. 

Here we write the vector as the superposition (with coefficients x, y, z) of three basis 
vectors. It means that the irrep F is contained thrice in the decomposition of the mechanical 
representation into irreducible representations of the group G = [a^]. 

This result can be generalized to the case of arbitrary and m in trivial manner: each 
irrep of the group G = [a™] enters exactly m times into the decomposition of the mechanical 
representation for A^-particle chain, and the rule for constructing m appropriate basis vectors 
is fully obvious from Eq. 

C. Proof of TheoremlH 

Proof. The basis vectors of all irreps Fj, listed for the case = 12, m = 3 in Table|lJ can be 
obtained from (jHHjl setting 7 = l,i, —1, —i, respectively (these values are one-dimensional 
matrices corresponding in F, (i = 1, 2, 3, 4) to the generator g = a^). 

Let us write the above basis vectors sequentially, as it is done in Table lTTl and prove that 
12 X 12 matrix, determined by this table, is precisely the matrix S that splits the linearized 
dynamical equations d = J{t) ■ d for the considered case. In TablelTTl we denote the basis 
vectors 0j(Fj) by the symbol of the irrep Fj [i = 1, 2, 3, 4) and the number j = 1, 2, 3 of the 
basis vector of this irrep. The normalization factor (|) must be associated with each row of 
this table to produce the normalized basis vectors (because of this fact, we mark the rows 
as 2 ■ (pj{Ti) in the last column of Table ITT|l . 

Obviously, we can use the matrix S from TablelTTl (the rows of this matrix are the basis 
vectors of all the irreps of the group T4) not only for the action on the vectors in the X- 
space of the full nonlinear system, but on the vectors in the (5-space of the linearized system 
6 = J(t) ■ 6, as well. It is essential that in the latter case the matrix S reduces the Jacobi 
matrix J(t) to a certain block-diagonal form. Indeed, as was shown in Theorem^ the matrix 
J(t) commutes with all the matrices of the mechanical representation of the bush symme- 



try group. Therefore, according to Wigner theorem 
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17|, it can be reduced, using unitary 



TABLE II: Basis vectors of the irreducible representations of the cychc group T4 
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transformation by the matrix S, to the block-diagonal form with blocks whose dimension is 
equal to rij ■ rrij. Here rij is the dimension of the irrep Tj, while rrij is the number of times 
that this irrep enters into the decomposition of the mechanical representation (constructed, 
in our case, in the (5-space). 

In Sec. lIVBl we have shown that for the translational group T = [a™] all rij = 1 and all 
rrij = m. Therefore, the above matrix S decomposes the Jacobi matrix J(t) into blocks whose 
dimension is equal to m. As a consequence of this decomposition, the system S = J(t) ■ 5 
splits into k = N/m independent subsystems, Lj {j = 1,2,..., A;), each consisting of m 
differential equations of the second order. The coefficients of these equations are time- 
dependent functions, and this time dependence is determined by the functions A{t), B{t), 
C(t) etc., entering into the bush displacement pattern (see, for example, ()7f8|l ). For the 
one-dimensional bushes, the coefficients of the above subsystems Lj turn out to be periodic 
functions with identical period, while for the many- dimensional bushes they possess different 
periods (such bushes describe quasiperiodic motion). 
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In general, it is impossible to obtain the explicit form of the functions A{t), B{t), C{t) 
etc., determining the bush displacement pattern. Therefore, we can add the bush dynamical 
equations to the m differential equations of each subsystems Lj. These additional d equations 
determine the functions A{t), B{t), C{t) etc. implicitly, where d is the dimension of the 
considered bush B[a™, . . . ]. 

On the other hand, we can give the following estimate for the bush dimension d: 

d<m. (59) 

Here m is the index of the translational symmetry of the bush B[d™, . . . ] (it determines the 
ratio between the size of the primitive cell in the vibrational state and in the equilibrium). 
Indeed, the bush displacement pattern can be found as the solution to the equation GX = 
X. If we take into account only translational symmetry group of the bush B[a'", . . .], i.e. 
G = [a"^], this equation reduces to the equation g(f) = cf) {g = d™) for the basis vector (f) of 
the identity irrep (7 = 1 in ((SEI)) of the group T = [d™]. As it has been already shown in 
Sec. lIVBl such vector (f) depends on exactly m arbitrary parameters. But some additional 
symmetry elements, denoted by dots in the bush symbol B[d™', . . .], can lead to a decrease 
in the number m of the above parameters. For example, the bush B[d^,i] turns out to be 
one-dimensional, i.e., in this case, d = 1, while m = 3. Even the vibrational bush B[d^] 
turns out to be two-dimensional (d = 2), if the condition of immobility of the mass center 
is taken into account. Thus, for all cases, d < m, and we can state that m equations of 
each Lj, extended by d additional equations of the given bush, provide us with k = N/m 
independent subsystems Lj of (m + d) autonomous differential equations. 

Taking into account the additional bush symmetry elements, denoted by dots in the 
symbol B[d'", . . . ], leads not only to reducing the bush dimension, but to a further splitting 
of the above discussed subsystems Lj (we consider this point in the next section). Thus, 
the linear stability analysis of the bush B[d'", . . . ] in the A^-particle chain indeed reduces to 
studying stability of individual subsystems whose dimension does not exceed {m + d). This 
is the conclusion of Theorem |21 and, thus, we have completed the proof. □ 
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D. Example 1: splitting the linearized system for the bush B[a^] 



We consider the sphtting of the hnearized system 5 = J{t) ■ 6 for the bush B[a^] in a 
chain with = 12 particles. The original nonlinear system, for this case, reads 



i - 



f{x. 



i+l 



1,2,..., 12 (xo = Xi2, Xi3 = Xi). 



(60) 



The displacement pattern of the bush B[a^], obtained from the equation d^X = X reads: 

X = {x{t),y{t), z{t) I x{t), y{t), z{t) \ x{t),y{t), z{t) \ x{t), y{t), z{t)}. (61) 

Substituting this form of vibrational pattern into ()6U|) . we obtain three differential equations 
for the functions x{t), y{t), z(t) (all the other equations of (plj) turn out to be equivalent to 
these equations): 

X = fiy - x) - f{x - z), 

y = f{z-y)-f(y-x), (62) 

z = f{x - z) - f{z - y). 
The linearization of the Eqs. (|6U|) near the dynamical regime determined by (j61|) leads to 
the system 

5 = J(t) ■ 6 (63) 

with the following Jacobi matrix: 
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(64) 
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where 



A{t) = f'[y{t) - 
B{t) = f[x{t)-z{t)], 

C{t) = f'[z{t)-y{t)l 

a{t) = -[Ait) + Bit)], 

f3{t) = -[A{t) + C{t)], 

j{t) = -[Bit) + C{t)]. 

Using TablelTTj the matrix S that sphts up the system (jH!^ can be written as follows: 



(65) 





/l 


I I 


A 


1 


I 


-il -I 


il 


2 


I 


-I I 


-I 






il -I 


-^l) 



(66) 



where I is the 3x3 identity matrix 



I 



1 
1 
1 



With the aid of the unitary transformation 

Jneiy(^) = • J(t) ■ S, 



(67) 



we obtain 



At) 



where 



V 



with 7i = 1, 72 = i, 73 



iA + B) 
A 
IkB 

■ -1, 74 ■ 



/ Di \ 
D2 
D3 

V D4/ 

A -fkB ^ 

-(A + C) C 

C -{B + C) j 
- —i ilk is the complex conjugate value of 7^). 



(68) 



(A; = 1,2,3,4), 



(69) 



This means that the hnear transformation 



(70) 
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reduces the old equations (jU^ to the following form 

with block-diagonal matrix Jnew{t) = S+ ■ J(t) • S determined by Eqs. ()68I69|) . 
Assuming 

we can present fl71|) in a more explicit form: 

5f ) = -(A + B)5f^ + A5f^ + iuB5f\ 



5^'=^ = A5r -{A+ c)5r + (72) 

where 7^ = 1, i, —1, — i for k = 1^2, 3, 4, respectively. 

Thus, we obtain four independent 3x3 systems of linear differential equations with 
time-dependent coefficients v4(t), B{t) and C(t), which are determined by Eqs. (jHSj). 

Let us write down these equations for the FPU-a chain. For this case, the function f{x) 
in (IHUjl reads f{x) = x + x^. Therefore, f\x) = 1 + 2x, and we obtain from (pTK|) 

A{t) = l + 2[y{t)-x{t)l 
Bit) = l + 2[x{t) - z{t)l 
C{t) = l + 2[z{t)~y{t)]. 

Substituting these functions into (f7^ one can finally obtain the following equations for the 
FPU-a chain: 

= hit) ■ 



where 



^ -2(1 + y~z) l + 2iy-x) jk[l + 2ix - z)]^ 



hit) 



(A; =1,2,3,4). (73) 



l + 2iy-x) -2il + z-x) l + 2iz-y) 
y%[l + 2ix-z)] [l + 2iz-y)] -2il + x-y) J 

Here 71 = 1, 72 = i, 73 = —1, 74 = —i and x = x(t), y = y(t), z = z(t) are functions 
determined by the dynamical equations of the bush B[a^] : 

X = iy — 2x + z) ■ il + y — z), 

y = iz -2y + x) ■ il + z - x), (74) 
z = ix — 2z + y) ■ il + X — y) . 
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These equations can be obtained from ()62j) taking into account the relation f{x) = x + 
for the case of the FPU-a modeL 

Remark. According to (|62j) (see, also ((73))), the vibrational bush B[a^] is three- 

dimensional. However, it actually turns out to be a two-dimensional bush. Indeed, there 
is no onsite potential in the FPU-a chains and, therefore, the conservation law of the total 
momentum of such system holds. Assuming that the center of masses is fixed, we obtain 
an additional relation x{t) + y{t) + z{t) =0 which reduces the dimension of the bush B[a^] 
from 3 to 2. 

E. Further decomposition of linearized systems based on higher symmetry groups 

Up to this point, we have discussed the decomposition of the linearized system 5 = J(t)-S 
using only the translational part of the bush symmetry group. In general, one can arrive at 
a more detailed splitting, if one takes into account the additional bush symmetries. 

1. Example 2: Splitting of the linearized system for the bush B[d^,i] 

Let us consider the decomposition of the linearized system for the bush B[a^,i] in the 
case of an arbitrary monoatomic chain. Since translational part of the symmetry group 
G = [a^,i], which turns out to be the dihedral group, is the same as that of the early 
considered bush B[a^], we can take advantage of all the results obtained in Sec. lIVE)1 and 
add only some restrictions originating from the presence of the additional generator i of the 
group [a^,i]. 

Substituting the vector X{t) in the form (jblj) into the equation iX{t) = X{t), we obtain 
z{t) = —x{t), y{t) = —y{t) and, therefore, y{t) = 0. The displacement pattern for the bush 
B[a^,i] then can be written as follows: 

X{t) = {x{t), 0, -x{t) I x{t), 0, -x{t) I x{t), 0, -x{t) I x{t), 0, -x{t)}. (75) 

Thus, the bush B[a^,i] turns out to be one- dimensional. 

As a result of the substitution z{t) = —xit), y{t) =0, three equations reduce to only 
one equation 

i = f{-x)- f(2x). (76) 
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For the FPU-a chain (see Eqs. ()74p ) . this equation transforms to 

X + 3x + = 0. 



(77) 



Unhke the purely translational group [d^], of the three-dimensional bush B[d'^], the sym- 
metry group [d^,i] of the one-dimensional bush B[a^,i] is the dihedral group with another 
sets of the irreps and basis vectors. It can be shown that taking into account that [d^,i] 
is the supergroup with respect to group [d^], allows one to obtain the following splitting 
scheme of the linearized system d = J{t) ■ S 

l:('5i); ^--m 2: (5,, 5,); ^^^^ 

2 : (5^,65); 3 : (67,69,511); 3 : {6s, 610, 612). 

Here we present the dimension of each independent subsystem (before the colon) and the list 
of its variables (after the colon). From the scheme ()78|). one can see that two of four three- 
dimensional subsystems corresponding to the splitting provided by group [a^] (see Eqs. fl7^ ) 
in the case of the supergroup [d^,i] are decomposed into new independent subsystems of 
dimensions equal to 1 and 2. Below we explain how one can obtain the splitting schemes 
analogous to (ffHjl . 



2. General case 

Now we consider the application of the Wigner theorem in case of an arbitrary bush B [G] . 
Let us consider a matrix H commuting with all matrices of a reducible representation F 
of the group G that can be decomposed into the irreps F^ of this group as follows 

F = $:®m,F,. (79) 
j 

According to the Wigner theorem, the matrix H can be reduced to a block-diagonal form 
with the blocks of dimensions rij ■ rrij, corresponding to the each Tj, with rij being the 
dimension of the irrep Tj entering rrij times into the decomposition (f?^ of the representation 
F into the irreducible parts. 

Moreover, each block possesses a very specific form, namely, it consists of subblocks 
representing matrices proportional to the identity matrix of the dimension Uj which 
repeat rnj times along the rows and columns of the block Dj. We can illustrate the structure 
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of a certain block Dj = D characterized by the numbers nj = n, rrij = m as follows 



D 



/^2lln ^J'22^n ■ ■ ■ l^2mXn 



(80) 



where I„ is the n x n identity matrix. 

In our case, the matrix H is the Jacobi matrix J(t) of the linearized system S = J(t) ■ 6, 
G is the symmetry group of a given bush B[G'], F is the mechanical representation of this 
group. Each block D^- generates an independent subsystem with rij ■ ruj equations in the 
decomposition of the linearized system. However, each of these subsystems automatically 
splits into rij new subsystems consisting of rrij differential equations, as a consequence of the 
specific structure of the block (see Eq. ()80|) ). Indeed, for example, if a certain D-block 
for the matrix J(t) possesses the form {rij = 3, nij = 2): 



A^2ll3 fJ'22^3 



it is easy to check that we obtain the following three independent pairs of the equations 
from the system S = J{t) ■ S: 

6i = findi + 121264, 
64 = 1^2161 + 122264-, 



62 = 1^1162 + 1^126^, 

4 = 1^2162 + /i22^5, 



(81) 



^3 = l^ii6-i + 12126%, 
6& = 122163 + 122266. 

Note that the dimension of each subsystem (jHTjl is equal to rnj = 2, while the total number 
of these subsystems is equal to rij = 3. 
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3. Irreducible representations and their basis vectors for the dihedral group 



Hereafter, for simplicity, we will discuss only chains with an even number (A^) of particles 
and illustrate the main ideas with the example = 12. 

The symmetry of an A^-particle (monoatomic) chain is completely described by the di- 
hedral group Go = Dn which can be written as the union of two cosets with respect to its 
translational subgroup T/v = {e, a, a^, . . . , a^~^}: 

Dn = Tn®Tn- i. (82) 

Here i is the inversion relative to the center of the chain. The group is a non-Abelian 
group, since some of its elements do not commute with each other (for example, ia = a~H). 
As a consequence, the number of classes of conjugate elements of this group is less than the 
total number (2A^) of its elements and some irreps Tj are not one-dimensional. The irreps of 
the dihedral group can be obtained by the well-known induction procedure from those 
of its subgroup T/y. It turns out that for Df^f with even N there are four one-dimensional 
irreps, while all the other (y — 1) irreps are two-dimensional. We discussed the construction 
of these irreps in where the following results were obtained. 

Every irrep can be determined by two matrices Mj{d), Mj(i) corresponding to its gener- 
ators a and i, where j is the number of t lis irrep. Four one-dimension irreps {j = 1, 2, 3, 4) 



are real and are determined by matrices 

Mj{d) = ±1, Mj{i) = ±1. (83) 
All other irreps are two-dimensional and are determined by matrices 



M,(a) , M,(i) 





with fij = e^, fLj = e-^ (j ^ 0, A^/2) [29|. 

Let us find the basis vectors of the irreducible representations of the dihedral group 
[a™,!] for the case m = 3 which corresponds to the bush B[a'",i]. Let and if) be the basis 
vectors of the two-dimensional invariant subspace corresponding to the irrep with the matrix 
( 1 ^\ 

^{9) = ! where ^ = a™ is the translational generator of the dihedral group. They 

can be obtained from the equations g4> = and gtl) = •ytl), respectively. For example, using 
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Eq. (fSHj) for the case m = 3, we find = {x, y, z \ 7~^x, ■j~^y, 'y'^z \ 7~^x, 7~^?/, 7~^-2 | . . . }, 
ip = {x,y,z I ^~^x,'y~^y,j~^z \ 7~^x, 7~^y, 7~^5 | . . . }. Here {x,y,z) and (x, 5) are 
arbitrary constants which these vectors depend on. 

Taking into account the presence of the matrix M(i) = in every two-dimensional 

irrep, one can state that 

10 = V, iip = 4>- (84) 

Because of these relations, there appear certain connections between the arbitrary constants 
{x, y, z) and (x, y, z). As a consequence, the basis vectors cf) and i/?, for each two-dimensional 
irrep of the group [a^,i], depends on only three arbitrary parameters: x, y and z. In turn, 
this means that each two-dimensional irrep enters exactly three times into the decomposition 
of the mechanical representation of the dihedral group [a^,i]. 

Unlike this, the one-dimensional irreps of the dihedral group [a^,i] are contained in the 
mechanical representation less than 3 times. Indeed, let us consider the basis vectors cf) and 
i/? of the one- dimensional irreps of the group [d'^] determined by the matrices M(a^) = 1 and 
M(a^) = —1, respectively, for the case m = 3, = 12: 

(f) = {x,y,z \ x,y,z \ x,y,z \ x, y, z}, 

'il) = {x,y,z\ - x,-y,-z \ x,y,z \ -x,-y,-z}. 

These vectors can be obtained from Eq. ()58|) by letting 7 = 1 and 7 = — 1. If the vector 
4> is not only the basis vector of the irrep of the group [a^], but also is the basis vector of 
a certain one-dimensional irrep of the dihedral group [d^,i], it must satisfy the equations 
10 = (/), for the irrep Fi (M(i) = 1) and i0 = —0, for the irrep r2 (M(i) = —1). We obtain 
z = —X, y = from the former equation, and z = x for the latter equation. Thus 

4>[Ti] = {x, 0, —X I X, 0, —X I X, 0, —X I X, 0, —x}, 
(j)[T2] = {x,y,x \ x,y,x \ x,y,x \ x,y,x}. 

In the same manner, we obtain the basis vectors V'irs] and i/^lTi] from the equations 
1-0 = 1/) and 10 = —0, respectively: 

^] = {x,y,x\ - x,-y,-x \ x,y,x \ -x,-y,-x}, 
0[r4] = {x, 0, —X I — X, 0, X I X, 0, —X I — X, 0, x}. 
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From the above results, we conclude that the irreps Fi and F4 are contained once, while the 
irreps F2 and F3 are contained twice in the decomposition of the mechanical representation 
for the considered chain. 

The generalization of these results to the case of the dihedral group [a™,i] with arbitrary 
m is trivial. 

The splitting scheme (fTHj) for the bush B[a^,i] for the monoatomic chain with = 12 
particles can be now explained as follows. There are five irreps (Fi, F2, F3, F4, F5) of the 
group [a^,i] = {e,a^,d^,a^ | i, ia^, ia^, id^} = D^^. As have just shown, the one-dimensional 
irreps Fi {ni = 1) and F4 (n4 = 1) are contained once (mi = 1, = 1) in the decomposition 
of the mechanical representation F for our chain. On the other hand, the one-dimensional 
irreps F2 (^2 = 1) and F3 {n^ = 1) are contained twice (m2 = 2, = 2) in F, while 
two-dimensional irrep F5 {n^ = 2) is contained thrice (ms = 3) in F. The twelve variables 
6j from Eq. ()78p are associated with the irreps of the group ^4 in the following manner: 

5i-ri(l); 

62 ^ F2(l), ^3 ^ F2(2); 

5, F3(l), 6, -> F3(2); (85) 

^6 r4(i); 

{67, 6s) ^ F5(l), {6,, 610) ^ r5(2), {6u, 612) ^ r5(3). 

Here, in parenthesis, we give the index of the copy of the irrep Fj (whose dimension is 
equal to n^) in the decomposition of the mechanical representation F. Note that the total 
number of such copies determines how many times the irrep Fj is contained in F. On the 
other hand, as we already know, rrii shows us the dimension of the subsystems Lj, while 
determines the total number of Lj with the same dimension associated with Fj. As a result, 
we obtain the splitting scheme (fTSjl . 

The above discussed decomposition of the full linearized system S = J(t) ■ d into inde- 
pendent subsystems Lj of small dimensions permits one to analyze efficiently the stability 
of a given bush in the monoatomic chain with arbitrary large number of particles (A^). 

Using this idea, the stability diagrams for all the one- dimensional bushes in both FPU-a 
and FPU-/3 chains were obtained in 9]. As an example, in Fig. Q we reproduce the stability 
diagram for the bush B[d^,i] for the FPU-a chain from that paper. In this diagram, each 
point {A, q) determines a certain value of the bush mode amplitude A and a certain value 
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1 




FIG. 1: Regions of stability (white color) of different modes of the FPU-a chain, interacting 
parametrically with the one-dimensional bush B [a?' , i] . 

of the wave number q = ^ that is associated with the index j of a fixed mode. The black 
points (A, q) correspond to the case where the mode j = becomes excited because of 
its parametric interaction with the mode of the bush B[a^,i]. The white color denotes the 
opposite case: the corresponding mode j, being zero at the initial instant, continues to be 
zero in spite of its interaction with the considered bush. Such a diagram allows one to study 
stability of one-dimensional bushes not only for finite A^, but also for the case N oo (some 
more details can be found in j^). 



V. SOME ADDITIONAL EXAMPLES 

A. Stability of the vr-mode in the FPU-a chain 

The symmetry group of the vr-mode X{t) = {A{t), -A{t), A{t), -A{t), A{t), -A{t), ...} 
in the FPU-a chain is [a^,i]. Using results of the previous sections, one can deduce for 

divisible by 4 that the linearized system in the vicinity of this mode splits into four 
individual equations and a number of two-dimensional systems of differential equations 2^ . 
The first pair of individual equations represents two independent harmonic oscillators, while 
the second pair represents two Mathieu equations. All the other systems represent pairs of 
coupled equations: 

5, + 4sin2 (f ) 5, = r^sin (^) s^_^ cos(2t), 
6n_. + 4cos2 ) 6n_^ = r/sin (^) 6j cos(2t). 
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where r/ = j = 1, 2, . . . , ^ — 1. In (jHUj) . A is the amphtude of the vr-mode, i.e. A = 
max|A(t)|. Eqs. fj86|) can be rewritten as follows: 

X + 4sin^ (I) X = 1] sm{q)y cos{2t) , ^^^^ 
y + 4: cos^ (I) ?/ = ?7 sin(g)x cos(2t), 

where g = ^ is the wave number, x{t) = Sj{t), y{t) = 5N__-{t). These results were obtained 
and discussed in jsj using a different (in comparison with the present paper) method. Indeed, 
there we obtained exact equations for the FPU-a chain in the modal space and only then 
linearized them near the vr-mode. 

As it was discussed in the system (jSTjl turns out to be rather remarkable: the vr- 
mode (or one-dimensional bush B[a^,i]) loses its stability simultaneously with respect to 
interaction with all the other modes. In other words, the threshold r/c for the loss of stability 
for the bush B[a^, i] turns out to be the same for all the values of g, i.e. for all the subsystems 



r/c = 2.42332 . . . (88) 

This property of the linearized system (|87jl was discussed in more detail in . The stability 
diagram for the bush B[a^,i] can be found in j^. 



B. Stability of the vr-mode in the FPU-/? chain 

This case is of a particular interest since the linearized system for the vr-mode in the 
FPU-/5 chain possesses a higher symmetry compared to that in the FPU-a chain. Indeed, 
for the FPU-/3 model, the interparticle potential is an even function, and some additional 
symmetry elements of the dynamical equations appear as a consequence of this fact. 

Let us introduce an operator u that changes signs of the displacements of all the particles: 

u{xi,X2, . . . , Xn} = {-xi, -X2, . . . , -xn}- (89) 

This operator leaves the FPU-/3 Hamiltonian unchanged because of its evenness. 

Therefore, the symmetry group of the FPU-/5 dynamical equations turns out to be a 
supergroup G'q with respect to the sjTumetry group Gq = [a, i] of the FPU-a chain: 

G"o = Go © Go ■ M (90) 
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(note that = e). 

As a result, one can classify bushes of modes in the FPU-/3 chain by subgroups of the 
group G'q rather than by subgroups of the Gq. The invariant manifolds (bushes of modes) in 



the FPU-/? model with respect to the group Gq were found by Rink in ll| . Some additional 
details of this problem were discussed in our paper 9] (in particular, the dynamical equations 
and stability of these bushes of modes). 

Considering G"q from ()90|) as the parent group, we discover that the 7r-mode in the FPU-/? 
chain must be characterized by the group G = [d^, i, au], unlike the group [d^, i] characterizing 
the vr-mode in the FPU-a chain with the parent group [d, i]. 

The third generator |31|J ait = ua of the group G = [d^,i, d-d] acts on the configuration 
vector X as follows: 

d'dX = {-xat, -xi, -X2, . . . , -xn-i}- (91) 

Above, we specified the group G by three generators, but it can be determined by only 
two generators. Indeed, if 

p = au, (92) 

then = d-dd-d = d^-d^ = d^. Thus, the first generator in the list G = [d^,i, d-d] is simply 
the square of the third generator p = au. 

For simplicity, let us consider the case = 6. Then [d^] = {e, d^, d^}, (d^ = e) is a cyclic 
group of the order 3 and the full group G contains 3 • 2 • 2 = 12 elements. With the aid of 
the operators p = au and g = i, we can obtain all the elements of the group 

G = [d^,i, d-d] = [p,q] (93) 

as follows: 

G = n®Te-q, (94) 
where Tg is the cyclic group of the order 6: 

n = {e,p,p',p\p\p'}. (95) 

The following generating relations, fully determining the group G = [p,q], can be ob- 
tained: 

p^ = e, = e, qp = ff'q. (96) 
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TABLE III: Irreducible representations of the dihedral group Dq = [p,q]. In this table // = e 3 , 
u = fi = e~ , p, and u are conjugate complex numbers with respect to fx and v, respectively. 

M{p) M{q) 
1 1 
1 -1 
-1 1 
-1 -1 

1 

1 

1 

1 



From these relations, one can see that G is the dihedral group Dq (by the way, it is isomorphic 
to the point groups Cgt,, -D3/1 and D^^, as well). 

The irreducible representations of the dihedral group were discussed in (see also 
Sec. lIVEl of the present paper). There are four one- dimensional and (y — l) two-dimensional 
irreps of the group with an even index A^. For simplicity, let us discuss the case = 6 
only (the generalization to arbitrary values of N turns out to be trivial). 

All irreps of the group Dq can be constructed from the irreps of its subgroup Tq (see 
(I94p ) with the aid of the induction procedure. As a result, we obtain the following irreps of 
the group G = [p, q] presented in Table lTTTl In this table, each irrep Tj (j = 1,2,..., 6) is 
determined by the matrices M{p) and M(g), corresponding to the generators p and q of the 
group Dq. 

All the invariant subspaces of the configuration space of the chain with A^ = 6 particles, 
corresponding to these irreps, and their basis vectors can be obtained as follows. Let us 
find the basis vector of a certain one- dimensional irrep F determined by the matrices 
M{p) = 7 and M(g) = 6 from Table lTTTl (7 = ±1, 6 = ±1). Thus, the vector (p must satisfy 
the relations: 

pcf) = -fcj), (97a) 
qcf) = 6(1). (97b) 
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From ()97ap . we obtain 

= {x, —'j^x, 7^x, — 7^x, 7^x, — 7x}, (98) 
where x is an arbitrary constant. For 7 = 1, the basis vector from (|98|) transforms to 

4>i = {x, —x,x, —x,x, —x}, (99) 
while for 7 = — 1 it transforms to 

02 = {x, X, X, X, X, x}. (100) 

Now we must also demand ()97b|) to hold. Remembering that ^ = i, and, therefore, that 
q(f) = {—Xq, —X5, —X4, — X3, —X2, —Xi}, we obtain 

q4>i = {x, —X, X, —X, X, —x} = (1)01, 

q4>2 = {-X, -X, -X, -X, -X, -x} = (-1)02- 
This means that 0i from ()99p is the basis vector of the irrep Fi (7 = 1,5 = 1, see Eqs. (jHZI)), 
while 02 from ()100|) is the basis vector of the irrep F4(7 = — 1, 5 = — 1). Thus, each irrep 
Fi and F4 is contained once in the decomposition of the mechanical representations. 

On the other hand, the basis vectors of irreps F2 (7 = 1, 5 = — 1) and F3 (7 = —1, 
6 = 1) are equal to zero. Indeed, demanding g0i = (— l)0i (5 = —1) one concludes that 
X = —X and, therefore, a; = and 0i = 0. The same result originates also from the equation 
g02 = (1)02 (5 = 1), namely, 02 = 0. Thus, we must conclude that the irreps F2 and 
F3 are not contained in the decomposition of the mechanical representation of our chain 
into irreducible representations (the corresponding invariant subspaces turn out to be null 
spaces) . 

There are two basis vectors, and "0, for each two-dimensional irreps from Table lTTTl For 
the irrep F5 the following relations must hold 

p(f) = fi<p, (101a) 

pip = fiil^, (101b) 

g0 = V, (101c) 

qip = cf). (lOld) 

Comparing ()101a|) and ()101b|) with ()97a|) . we can write 

= {x, —fi^x, fi'^x, —fi^x, fi'^x, —fj,x}, (102) 
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V' = {y, -fj'^y, fj'^y, -f^^y, /^^y, -f^y}, (103) 

where x and y are two different constants. 

Then, from (jlOlcj] or (jlOldjl . we obtain a certain relations between the constants x and 

y- 

y = fix. (104) 

Substituting this value of y into ()103|) we, finally, conclude that the irrep is contained only 
once in the decomposition of the mechanical representation of our FPU-/? chain. Indeed, 
the basis vectors <p, if) turn out to depend on only one constant x which can be determined 
from the normalization condition. The same conclusion is valid also for the second irrep Fg 
of the group Dq. 

The generalization of this conclusion to the case of the FPU-/5 chain with an arbitrary 
even number = 2n of particles can be achieved trivially. Namely, each two-dimensional 
irrep of the group is contained only once in the mechanical representation of this group 
for the FPU-/3 chain. One-dimensional irreps are contained in the mechanical representation 
once or not at all. Remembering the above discussion about the application of the Wigner 
theorem to splitting the linearized system 5 = J{t) ■ 5, we discover that for the FPU-/5 chain 
(unlike FPU-a chain) this system is decomposed into individual differential equations of the 
second order, i.e. we have the complete splitting in this case. 

The following remark should be done to avoid a possible misunderstanding. The form of 
the TT-mode 

X{t) = {Ait), -Ait), Ait), -Ait), Ait), -Ait)}, (105) 

is one and the same 1^32] for both FPU-a and FPU-/5 chains. Therefore, its symmetry group 
can be written as [a^,i, a-u] rather than [a^,i] since the operator aii does not change the 
pattern p05|) not only for the FPU-/3 chain, but also for the FPU-a chain. 

Nevertheless, the linearized system 6 = Jit) ■ S, for the FPU-/5 chain, is invariant with 
respect to the operator au, while that for the FPU-a chain is not invariant under the action 
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of du. Indeed, the hnearized (near the 7r-mode) system for the FPU-a chain reads: 



Si = 


[62 


-26i 




- m) 


[S2 


-s. 


§2 = 


[S3 


-252 




+ 4:A{t) 


[63 




53 = 




-263 


+ 62] 


- m) 


[54 


-62 


s,= 




-264 


+ S3] 


+ 4A(t) 




-S3 


k = 


[^6 


-26, 


+ 54] 


- 4A(t) 




-h 


k = 


[^1 


-26e 


+ 55] 


+ 4:A{t) 


[5i 


-h 



(106) 



On the other hand, for the FPU-/3 chain, the hnearized system reads: 



51 = 


[S2 


- 26i 




■ [1 + 12A\t) 


82 = 


[S3 


-252 


+ 5i] 


■ [1 + 12 A\t) 


53 = 


[54 


-263 


+ S2] 


■ [1 + 12A\t) 


5,= 


[^5 


-264 


+ S3] 


■ [1 + 12^2 (t) 


k = 


[^6 


-26, 


+ 54] 


■ [1 + 12 A\t) 


k = 


[5l 


-26e 


+ 55] 


■ [1 + 12^2 (t) 



(107) 



The operator du acts on the vector 6 follows: 

dud = du{6i, 62, S3, 64, 65, 6q} = {-6q, -61, -62, -S3, -64, -65}. (108) 
The substitution 

6^ —Sr, §2 — ^ —S-\, S3 — » —S2, 

' ' (109) 

S4 -^3, ^5 -64, 6q — > -65, 

transforms, for example, the first equation of the system (jlOfij) to the form 

6e = [Si-26e + S,]-AA{t)-[6i-6,], 

which differs by the sign in front of 4:A{t) from the sixth equation of the system p06p . Con- 
trariwise, the transformation ()109|) leaves the system ()107p for the FPU-/? chain unchanged 
(it leads only to some transpositions of the individual equation in ()107|) ). 

What does it mean? Each bush (the vr-mode, in our case) is associated with a certain 
subgroup of the parent group Go, the symmetry group of the original nonlinear dynamical 
equations of the considered mechanical system. The operator u does not belong to the 
symmetry group Gq of the FPU-a chain (unlike the case of the FPU-/5 chain!) and, therefore, 
there are no subgroups of the group Gq whose elements contain u. Precisely this fact can 
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explain why we must not take into consideration the operator au for the case of the FPU-a 
chain, even through this operator does not change the pattern ()105|) of the vr-mode. 

The function A{t) from the displacement pattern (jl(J5j) of the vr-mode is determined by the 
dynamical equation of the one-dimensional bush B[a^,i]. As it was already discussed, this 
equation can be obtained by substituting the vector X{t) from Eq. (|1U5|) into the nonlinear 
dynamical equations of the FPU chain. For the FPU-a model this equation turns out to 
be the equation of a harmonic oscillator A + AA = 0, while for the FPU-/5 model it is the 
Duffing equation A + AA + IQA^ = 0. 

Thus, in both cases, the exact expression for the function A{t) can be found. For studying 
the stability of the vr-mode, we must substitute the corresponding expression for A(t) into 
the linearized system 5 = J{t) ■ 6. As we saw in the present section, the linearized system 
splits into individual equations for the FPU-/3 chain, while for the FPU-a chain it can be 
decomposed into four individual equations and (A^ — 4)/2 pairs of differential equations. All 
these equations turn out to be equations of the second order with time-periodic coefficients 
depending on the function A{t) jal^- 

Let us note that stability of the vr-mode in the FPU-a and FPU-/? chains was investigated 
by different methods in a large number of papers (see, for example. 
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U, but, to our best understanding, the influence of symmetry of these mechanical 
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models on stability analysis was not discussed. Unlike the above cited works, the bush 
stability analysis presented in this paper based on the symmetry-related arguments only. 
Therefore, our conclusion about the difference in splitting scheme of linearized systems for 
the vr-mode in the FPU-a and FPU-/5 chains can be automatically extended to all the other 
nonlinear chains with the same symmetry characteristics. In particular, we can conclude 
that the splitting of the linearized system into individual equations can be performed not 
only for the FPU-/3 chain, but for every chain with an even potential of the interparticle 
interaction. In contrast, it is impossible for the chains with arbitrary potential, and the 
FPU-a model is a simple illustration of this proposition. 

We would like to focus on the paper 3|, where some analytical results were obtained 
for the stability of the vr-mode in the nonlinear chains with a general even form of the 
interparticle interaction potential. The author of ^| succeeded in his analysis thanks to the 
decomposition of the linearized system into individual equations, but such analysis cannot 
be extended to the FPU-a chain precisely because the potential of this model is odd. 
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TABLE IV: Irreducible representations of the dihedral group D^. Here = e~ , Ji = e~~ = a , 



I. 



Real form 




C. Stability of two-dimensional bush B[a^,i] in the FPU-a chain 

Now, let us consider the stability of the two-dimensional bush B[d^,i] that can be deter- 
mined by the displacement pattern jsj 

X{t) = {A{t),B{t),-B{t),-A{t) I Ait), Bit), -B{t), -Ait) \ -A(t)} 

(110) 

(for simplicity, we start with the case = 12). 

The symmetry group of this bush is the dihedral group with translational subgroup 

T3 = {e,a^a®} (a^' = e). (Ill) 

This non-Abelian group (id^ = d^i) consists of six elements determined by the equation 

D3=T3©T3-i, (112) 

and possesses the following three irreps presented in Table |TV1 (we give there the two- 
dimensional irrep not only in the complex form, but in the real form also). 

The basis vectors of the one- dimensional irreps of the group T3 (|111|) can be found from 
the relation 

pcf> = 70, (113) 
where 7 = 1, /x, Ji. From this equation we obtain 

= {Xg, Xio, Xii, X12 I Xi, X2, X3, X4 I X5, X6, X7, Xg} = 

7{Xi, X2, X3, X4 I X5, Xg, X7, Xg I Xg, Xio, ^n, X12} (114) 
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and then 

4> = {Xl, X2, X3, X4 I "yXi, 7X2, 7X3, 7X4 I 7X1, 7X2, 72^3, l^^}- (115) 

This basis vector depends on four arbitrary constants (xi, X2, X3, X4) and, therefore, it 
determines a four- dimensional subspace invariant under the translational group T3 (jlllj) . 
associated with the irrep F defined by the one-dimensional matrix M(p) = 7. 

The basis vectors of one- dimensional irreps of the whole group -D3 from pi2j) can be 
obtained with the aid of the equation icf) = (f) (for the irrep Fi) and with the aid of the 
equation Icf) = —(p (for the irrep F2), where 4> = |xi , X2, X3, X4 | xi, X2, X3, X4 | xi, X2, X3, X4} 



since 7 = 1 for these both irreps (see Eq. (jll5|) ) 



33| . From the equation lep = (p, we obtain 



0[Fi] = {Xi, X2, -X2, -Xi I Xi, X2, -X2, "Xi | Xi, X2, -X2, -Xi}. (116) 

This result means that the irrep Fi is contained twice in the decomposition of the mechanical 
representation of the considered chain. Analogously, for the case of the irrep F2, we obtain 
from the equation i0 = —0: 

0[F2] = {xi,X2,X2,Xi I Xi,X2,X2,Xi | Xi , X2 , X2 , Xi } . (117) 

This vector also determines the two-dimensional subspace and, therefore, the irrep F2 of the 
group is contained twice in the decomposition of the mechanical representation of our 
chain. 

Now let us consider the basis vectors of the two-dimensional irrep F3 from Table HVl The 
following relations must hold (for the complex form of this irrep) 

pcj) = (118a) 

pip = fill), (118b) 

l(p = -0, (118c) 

lip = 0, (118d) 

where (p and if) are two basis vectors of the irrep F3. Comparing ()118aj) and ()118bj) with 
(|115|) and letting 7 = /i and 7 = /i, respectively, we obtain: 

(p = {xi,X2,X3,X4 I pXi, flX2, pXs, pX4 \ p^Xi, p'^X2, p'^X^, p'^X^} (119) 

and 

V' = {Vh 2/2, 2/3, 2/4 I /^2/i, /^2/2, /^2/3, m I yU^yi, /^^2/2, /^^2/3, /^V}- (120) 
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On the other hand, the both equations ()118cp and ()118d|) lead to the same relation 
1(f) = if). Then, from Eqs. ()119|) . ()120p . we obtain yi = —^1x4, 1/2 = —fJ^x^, = ~fiX2, 
1/4 = — /iXi. Therefore, the final forms of the basis vectors of the irrep are: 

(f)[T-i\ = {xi,X2,X3,X4 I fLXi, fiX2, fiXs, J1X4 \ Jx^ Xi, Jl^ X2-, J^^ X^-, Ji^ X4} , (121) 



l/'lFs] = {-/ia;4, -/iX3, -/iX2, -/iXi I - f/x4, -/M^Xs, -fl'^X2, -fi^Xi, I - X4,-X3,-X2,-Xi}. 

(122) 

Each of these vectors depends on 4 arbitrary parameters (xi, X2, X3, X4) and, therefore, one 
can construct four independent pairs of the basis vectors — 0j[r3], ■j/'iirs] (j = 1,2,3,4) 
— of the irrep js^. In turn, this means that the irrep Fs is contained four times in the 
decomposition of the mechanical representation of the chain with = 12 particles. 

Taking into account the above results (see ()116p . (jllTj) . (jl2H) . (jl22p ). we can conclude 
that in the case = 12, the linearized system S = J(t) ■ S for studying the stability of the 
bush B[a'^, i] splits into two 2x2 and two 4x4 independent systems of differential equations 
of the second order. In the real form, these systems can be written as follows: 

i)i + K3{t)u, = Ki{t)u2, 
i)2 + K4{t)u2 = Ki{t)ui, 

Us + Ki{t)u3 = Ki{t)iy4, 
i)4 + K,{t)u4 = Ki{t)us. 



i>5 + KQ{t)u5 

1)7 + K^{t)iJj 
^ i>8 + Ki{t)us 



Ki{t)ue + K2it)ur, 



(123) 



^9 + KG{t)ug = Ki{t)uio + K2{t)uu, 

i)lO + K-i{t)uiQ = Ki(t)ug, 

i)u + K^{t)Vii = Ki{t)ui2 + K2{t)ug, 

^ i)i2 + Ki{t)ui2 = Ki{t)uu, 
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FIG. 2: Stability diagram for the bush B[a'*, i] in the FPU-a chain. 



where 

Ki{t) = 1 - 2A{t) + 2B{t), 
K2it) = V3[-l + 2A{t)], 
K^{t) = 3 + QA{t) + 2B{t), 
Ki{t) = 3-2A(t) -65(t), 
K,{t) = \-2A{t)-AB{t), 
K,{t) = l-2A{t). 

The generahzation to the case of an arbitrary (note that mod 4 = must hold!) is 
trivial: each one-dimensional irrep is contained twice in the decomposition of the mechanical 
representation of the A^-particle chain, while each two-dimensional irrep is contained 
four times in it. 

In conclusion, in Fig.0we reproduce the stability diagram for the two-dimensional bush 
B[a^,i] from our paper ||9|. This diagram corresponds to the FPU-a chain with = 12 
particles. It represents a planar section of the four- dimensional stability domain in the 
space of the initial conditions z/i(0), z/2(0), z>i(0), z>2(0), where viif) and z/2(t) are two modes 
of the considered bush. 
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In producing Fig.|21 we specify z^i(O) = 0, z>2(0) = 0, and change i^i(O), z/2(0) in some 
interval near their zero values. The stability domain, resembling a beetle, is drown in black 
color in the plain z/i(0), z/2(0). The bush B[a^,i] losses its stability (and transforms into 
another bush of higher dimension), when we cross the boundary of the black region in any 
direction. From Fig.|21it is obvious, how nontrivial stability domain for a bush of modes can 
be. 

A detailed description of the stability domains for one- dimensional and two-dimensional 
bushes of modes in both FPU-a and FPU-/? chains can be found in ^. 

VI. CONCLUSION 

All the exact dynamical regimes in A^-particle mechanical system with discrete symmetry 
can be classified by the subgroups Gj of the parent group Gq, i.e. the symmetry group of its 
equations of motion. Actually, each subgroup Gj singles out a certain invariant manifold 
which, being decomposed into the basis vectors of the irreducible representations of the 
group Go, is termed "bush of modes" 

The bush B[Gj], representing an ra- dimensional vibrational regime, can be considered 
as a dynamical object characterized by its displacement pattern of all the particles from 
their equilibrium positions, by the appropriate dynamical equations and the domain of 
the stability. One-dimensional bushes are symmetry-determined similar nonlinear normal 
modes introduced by Rosenberg jl6| (see also Q]). For Hamiltonian systems, the energy of 
the initial excitation turns out to be "trapped" in the bush and this is a phenomenon of 
energy localization in the modal space. 
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7|. Bushes of 
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3. 



The different aspects of the bush theory were developed in 
vibrational modes (invariant manifolds) in the FPU chains were discussed in 

The stability analysis of a given bush B[G'] reduces to studying the linearized (in the 
vicinity of the bush) dynamical equations S = J{t) ■ 5. In the present paper, we prove 
(TheoremHj) that the symmetry group of the linearized system S = J(t) • S turns out to 
be precisely the symmetry group G of the considered bush B[G]. This result allows one to 
apply the well-known Wigner theorem about the specific structure of the matrix (J(t), in our 
case) commuting with all the matrices of a fixed representation (mechanical representation, 
in our case) of a given group. According to the above theorem one can split effectively 
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the linearized system 5 = J(t) ■ d into a number of independent subsystems of differential 
equations with time-dependent coefficients. 

We want to emphasize that this symmetry-related method for splitting the linearized sys- 
tems arising in the linear stability analysis of the dynamical regimes is suitable for arbitrary 
nonlinear mechanical systems with discrete symmetry. Such a decomposition (splitting) of 
the linearized system S = J{t) ■ 6 is especially important for the multidimensional bushes of 
modes, describing quasiperiodic vibrational regimes, which cannot be treated with the aid 
of the Floquet method. Indeed, in this case, we need to integrate the differential equations 
with time-dependent coefficients over large time interval, unlike the case of periodic regimes 
where we can solve the appropriate differential equations over only one period to construct 
the monodromy matrix. 

The above method is applied for studying the stability of some dynamical regimes (bushes 
of modes) in the monoatomic chains. For this specific mechanical systems, we prove Theo- 
rem|21which allows one to find very simply the upper bound of dimensions of the independent 
subsystems obtained after splitting the linearized system 6 = Jit) ■ 6. Indeed, according to 
this theorem, the dimension of each such subsystem does not exceed the integer m deter- 
mining the ratio of the volumes of the primitive cell of the chain in the vibrational state, 
corresponding to the given bush B[G] =B[a"*, . . .], and the equilibrium state. 

Taking into account any other symmetry elements of the considered bush allows to reduce 
the dimensions of at least some of the above discussed subsystems. We illustrate this fact 
comparing the stability analysis of the vr-mode (zone boundary mode) for the FPU-a and 
FPU-/5 chains. 
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